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We investigate the nearest level spacing statistics of open chaotic wave systems. To this end we 
derive the spacing distributions for the 3 Wigner ensembles in the one-channel case. The theoretical 
results give a clear physical meaning of the modifications on the spacing distributions produced by 
the coupling to the environment. Based on the analytical expressions obtained, we then propose 
general expressions of the spacing distributions for any number of channels, valid from weak to strong 
coupling. The latter expressions contain one fit parameter. The surmise is successfully compared 
with numerical simulations of non-Hermitian random matrices and with experimental data obtained 
with a lossy electromagnetic chaotic cavity. 
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It is by now well established that classical chaos man- 
ifests itself in generic statistical properties of the corre- 
sponding quantum or wave systems. In the ideal case of 
closed systems, i.e. those whose coupling to the environ- 
ment can be neglected, the spectral and spatial statis- 
tics are well described by the Random Matrix Theory 
(RMT) pp. The statistics of chaotic wave systems co- 
incide with the Gaussian Orthogonal Ensemble (GOE) 
if time-reversal symmetry (TRS) holds; with the Gaus- 
sian Unitary Ensemble (GUE) if TRS is broken; and with 
the Gaussian Symplectic Ensemble (GSE) for TRS sys- 
tems with spin-1/2 interactions. Among all the statis- 
tical quantities used to analyze the spectral properties 
of closed systems, the nearest level spacing distribution 
P(s) is certainly the most referred one ( see Ref. [2] for a 
review). The spacings distributions, surmised by Wigner 
using the 2- level approximation [3], reads: 

P™ S ( S ) oc gPe^ , (1) 

where (3 is the Dyson index labeling GOE (/? = 1), GUE 
(/? = 2), and GSE (J3 = 4). These distributions were first 
used to describe spectral statistics of heavy nuclei and 
later successfully employed to describe spectral statistics 
of a wide range of closed and weakly open chaotic systems 
(see [1] and references therein). 

However, for systems of current interest such as quan- 
tum dots [5], nuclear compounds reactions 6J, micro- 
lasers cavities [7J [5] or microwave billiards [2] [TjHTT] , the 
coupling to the environment must be explicitly taken into 
account. The open systems arc then characterized by a 
set of resonances embedded in the continuum given by 
the poles of the S'-matrix [T2]. The poles are, in turn, 
the complex eigenvalues of the effective Hamiltonian: 

n eS = H- l -VV^ , (2) 

where the Hermitian part H is the Hamiltonian of the 
closed system giving rise to N real energy levels and 
the anti-Hermitian part iVV^ /2 models the coupling to 



the environment in terms of M scattering channels. The 
N x M matrix V contains the coupling amplitudes V£ 
that connect the nth level to the cth scattering chan- 
nel. The eigenvalues of the effective Hamiltonian are 
complex: £ n = E n — §r n , where E n and T n are, re- 
spectively, the eigenenergies and the resonance widths 
of the open system. Applying RMT to the effective 
Hamiltonian ^ (see [T3l [14] for recent reviews), the 
Hermitian part of % c ff is described by a Gaussian en- 
semble of the appropriate symmetry and the coupling 
amplitudes are considered as independent random Gaus- 
sian variables [15] : real for GOE, complex for GUE, and 
real quaternion for GSE, with zero mean and covariance: 

(VZ{V£,)*) = (l/r])6 nn ,5 cc ' . Here 1/r, is the coupling 

strength. 

Using the above formalism, much progress has been 
achieved in understanding the statistical properties of 
the widths from weak [16] to strong coupling [17j [18] as 
well as the statistics of eigenvectors [THl [2D] ■ However, to 
the best of our knowledge, there are still not general an- 
alytical expressions for the distributions of the spacings 
s„ = E n+ i — E n for open chaotic systems. This is so, even 
if analytical results for interesting particular cases have 
been reported [2"TH2"3"] and progresses have been done in 
analyzing the spacing distribution of the resonances in 
the complex plane [HJ [25] . What is missing then is the 
counterpart of the Wigner surmise for open chaotic wave 
systems. In this Letter, we achieve this goal by deriving 
analytically the probability distributions of the spacings 
for the 3 Gaussian ensembles for the 2-level model and 
the one-channel case. Then, we extend those results to 
the iV-lcvcl model and to any number of channels con- 
sidering the coupling strength as a free parameter. 

To derive the spacing distributions V^j =1 (s) we start 
with the joint energy distribution Pj^ =1 ({E n }, {r„}), 
first obtained for GOE by Sokolov and Zelevinsky [25] 
and then for the 3 ensembles by Stockmann and Seba 
[21] . In analogy with the derivation of the Wigner sur- 
mise, we assume that the 2-level approximation holds. 
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Specializing Eq. (4.4) of 
tion reads: 



the joint energy distribu- 
(E, E 2 ) 2 + Q±=pt 



(ri+r 2 )2 



{E 1 - E 2 y , 
.2=? _4^ 2 ^R 2 -£i£2)-^(r 1 +r 2 ) 



x {T 1 T 2 )^ A e- A ^ +E i-^-^ 



(3) 



where A fixes the mean level of the closed system: A = 
it/V^NA. X] A ^> 1 and rj A <C 1 respectively, correspond 
to the weak and strong coupling regimes. Changing to 
variables z — E% + E\ and s — E% — E\ in ([3| and inte- 
grating over z one gets: 



(ri-r 2 ) 



(ri+r 2 ) 2 
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x (r!r 2 )^ e - 



}(Ti+r 2 ) 



(4) 



To go further let us consider the 3 ensembles separately. 

For the GOE case, introducing in Q the new variables 
x = Ti + r 2 and y = I\ — T 2 , the integration over y can 
be done. With the proper normalization constant, the 
distribution yields: 



where /„ is the modified Bessel function of first kind. 

For the GUE case, the exponential prefactor can be 
written in terms of the partial derivatives of A and ?y, 
where TiT 2 and Ti + T 2 are generated by —2^ and 
respectively. The distribution can then be written 

as: 



^M= 2 i(^ r i> r 2)«e-¥- 
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drf dA 



x 6-^^-1(1*1+^) 



(6) 



The integrations over Ti and T 2 are now straightforward 
and give rise to: 



oo oo 



dTi / dT 2 e 



-4r 1 r,- m+ r,) = 2^ Ei (|_) j (?) 
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where Ei(a) = f°° dx^r is the exponential integral. 
Operating (J7J) with ^3 +2^- and including the nor- 
malization constant we finally arrive at: 



Ar\ _, 
16 e ~' 



dx 



(5) 



Z7T 



E(A,r)) s 2 



2 £(A )f? ) 
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(8) 

where 77) = Ei (^-J . Applying exactly the same 
method to GSE, one gets: 



A 3tt 2 3 
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E(A, V ) 



(9) 



Note that the distributions (5|, ^ and ^ tend to the 
Wigner surmise ([I]) in the limit of vanishing coupling 
1/rj — > 0. The Gaussian tail e _j4s / 2 , proper of the closed 
systems ([l]) , remains and the Wigner prefactor s* 3 now be- 
comes some function fp of s: V^ =1 (s) — fp(s)e~ As I 2 . 
These expressions reveal also that the effects of the cou- 
pling show up only in the prefactor. Thus, the main mod- 
ifications occur at small spacings. In particular, there ap- 
pears a finite probability for null spacings: V^ 1=1 {s) 7^ 0. 
In other words, the level repulsion, a main feature of the 
Wigner surmise, is suppressed for open systems. This ef- 
fect is associated to the phenomena of attraction of the 
levels along the real axis induced by the coupling to the 
environment [57] . For the GUE case fp=2 is a polyno- 



mial containing the Wigner term s 2 and a constant term. 
For the GSE case, in addition to the Wigner term s 4 and 
the constant term, a quadratic term appears too. The 
presence of this quadratic term can be viewed as a conse- 
quence of the breaking of the TRS, since it characterizes 
the behavior at small spacings for closed systems with 
broken TRS. 

It is important to note that in the limit of strong 
coupling I/77 — > 00 we obtain a Gaussian distribution 

^m=i( s ) — * \f^ e ~^ S ■ This behavior, intrinsic of the 
2-level approximation, does not correspond to the N- 
level model where the Wigner surmise is recovered |21j . 
Indeed, the 2-level model expressions are valid only up 
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to the intermediate coupling regime. However, we will 
show that if r/ is considered now as a free parameter, the 
distributions ([S]) and Q are valid for any coupling 
strength and, moreover, for any number of channels. 





FIG. 1: Distribution of the normalized spacings for the 3 
Gaussian ensembles with (r) /A = 1. On the top for GOE, 
on the middle for GUE, and on the bottom for GSE. On the 
left for M = 1 and on the right for M = 10. For the one- 
channel case, the exact analytical expressions pi), (JsJ> and Q 
are represented; for M = 10, considering 7/ as a fit parameter 
determined by a mean square fit procedure. In insert the tail 
of the distributions are represented in semi-log scale. For all 
curves, the mean level spacing is normalized to one. 



To confirm the validity of our proposal and to analyze 
how the parameter 77 evolves with the coupling and the 
number of channels, we carried out numerical simulations 
of non-Hermitian random matrices defined by j2|. The 
simulations were performed with ensembles of 150 non- 
Hermitian of size 1000 x 1000 varying the mean width 
(r) in the range [0.1, 30] A and for M = 1, 3, 5 and 10. 

A comparison between the expressions (|5|, d8k and ([§ 
and the numerical results is presented in Fig 1 For all 
coupling values and number of channels analyzed the con- 
fidence level is larger than 99.5%. We stress that the 
excellent agreement with the numerical simulations con- 
firm that the Wigner generalizations Eqs. ([5]), pi) and 
([9]), being exact for the 2-level model and M = 1, are ex- 
cellent effective formulas for arbitrary coupling strength, 
number of levels and channels. 

In Fig. [2]we show the evolution of the effective coupling 
strength 1 jr\ for the GOE case as a function of the mean 
width. We can distinguish two regimes. For weak enough 
coupling ((r) < 2A), 1/rj increases quite rapidly with the 
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FIG. 2: Evolution of the parameter r\ as a function of the 
normalized mean width (r) /A for the GOE case. M = 1 in 
solid line, M — 3 in dots, and M = 5 in dashed line. 



mean width. For strong coupling, the reverse occurs and 
we go back to the Wigner surmise as predicted in [21] for 
the one-channel case. Note, the approach to the Wigner 
surmise also occurs for M > 1 (with M <C N), albeit 
extremely slow. A similar behavior was observed for the 
GUE and GSE cases. 





FIG. 3: Distribution of the spacings for an open chaotic mi- 
crowave cavity. In histogram the experimental distribution, 
on the top with the range 1 to 6 GHz (M w = 10) and on 
the bottom 15 to 16 GHz (M w — 20). Both histograms have 
been obtained using 2600 events. The dotted line shows the 
Wigner surmise, and the solid line the best fit using expres- 
sion [5] The confidence level is lower than 50% for the Wigner 
surmise |l]) and larger than 90% for the expression |5|. 



Furthermore, we compare the theoretical prediction for 
GOE, Eq. ([5]), with recent experimental data obtained 
for an open chaotic 2D microwave cavity displaying TRS 



4 



in the closed limit. The data were provided to us kindly 
by U. Kuhl [18]. In that experiment, the coupling to the 
environment is produced by the finite conductivity of the 
wall (modeled by M w fictitious open channels) and one 
antenna attached to the cavity to perform the measure- 
ments. Fig. [3] shows the experimental distribution of 
the spacings for two different ranges of frequency i.e., for 
two different coupling strengths and number of coupling 
channels. The agreement is good with a confidence level 
larger that 90% in both cases. 

In summary, we have proposed general expressions for 
the distributions of the nearest level spacings of open 
chaotic wave systems, where the internal Hamiltonian 
belongs to one of the three well-known classes of Gaus- 
sian Ensembles. We stress again that the expressions are 
exact for N = 2 and M = 1 and otherwise excellent ap- 
proximations considering the coupling strength as a free 
parameter. The expressions indicated that coupling to 
the environment modified, with respect to the Wigner 
surmise, the form of the cofactors from s* 9 to some func- 



tion of s. Our predictions, valid from the weak to the 
strong coupling regime were shown to be in excellent 
agreement with numerical simulations of non-Hermitian 
random matrices and with experimental data from a mi- 
crowave cavity. In essence, the set of distributions |5|, 
^ and (|9| may be viewed as the extension of the Wigner 
surmise to open chaotic systems. 
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